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NYO-2883 
Abstract 

The  following  report  is  concerned  with  the  possible 
existence  of  a  shock  In  a  high -temperature,  colllsionless 
plasma  situated  in  a  m.agnetic  field.   The  underlying  equa- 
tions are  those  of  conservation  of  mass,  momentum  and  energy 
for  each  component  coupled  with  Maxwell's  equations.  The 
stress  tensor  is  assumed  to  be  anisotropic,  its  exact  form 
being  derived  from  guiding-center  theory.  The  conclusion 
Is  that  no  shock  solution  is  contained  in  this  system  of 
•equations. 
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NYO-2885 

ON  PLASMA-MAGNETIC  SHOCKS 

1.  Introduction. 

A  fully  ionized  fluid,  subjected  to  a  magnetic  field, 
can  be  expected  to  sustain  a  steady-state  shock  wave  if 
the  usual  damping  mechanisms  such  as  viscosity  and  heat 
conduction  are  present.  This  is  the  case,  for  example, 
at  low  temperatures  where  the  mean  free  path  for  collisions 
is  small  compared  to  the  dimensions  of  interest  and  hence 
the  concepts  of  viscosity  and  heat  conduction  are  meaning- 
ful. This  is  no  longer  the  case,  however,  at  temperatures 
sufficiently  high  so  that  the  mean  free  path  becomes  larger 
than  all  significant  dimensions.  Under  these  conditions, 
the  possible  existence  of  a  shock  is  not  at  all  obvious  and, 
in  fact  has  not  yet  been  proved. 

In  this  report,  we  shall  examine  the  properties  of  a 
set  of  equations  appropriate  for  a  high  temperature  plasma 
in  a  magnetic  field.  The  usual  Navier-Stokes  stress  terms 
which,  in  an  ordinary  fluid,  are  responsible  for  dissipation 
will  be  replaced  by  terms  derived  from  guiding  center  theory. 
We  shall  show  that  these  equations  do  not  possess  a  shock 
solution  and  that  the  only  possible  transition  is  a  pulse. 

Collision-dominated  or  "magnetohydrodynamic"  shocks  as 
they  are  usually  called  are  the  natural  starting  point  for 
investigating  shock  phenomena  in  ionized  fluids  since  they 
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more  closely  resemble  the  well-known  hydrodynamical  shocks. 
The  earliest  studies  were  concerned  with  deriving  the  jump 
conditions  and  led  co   the  discovery  of  three  possible  shock 
speeds.   Although  a  fully  ionized  fluid  is  composed  of 
positively  and  negatively  charged  particles  (ions  and 
electrons)  in  equal  numbers,  these  investigations  were  all 
based  on  a  one-fluid  approximation;  this  led  to  manageable 
equations  whose  mathematical  structure  closely  resembled 
that  of  ordinary  fluid  dynamics.   Moreover,  a  one-fluid 
theory  is  a  reasonable  approximation  at  sufficiently  low 
temperat\Ares  and  high  densities  to  enable  the  ions  and 
electrons  to  reach  thermodynamic  equilibrium  through 
collisions. 

The  problem  of  finding  a  smooth  transition  curve  joining 

the  constant  state  ahead  to  the  different  constant  state 

2 
behind  the  shock  was  first  considered  by  Marshall.   This 

again  was  a  one-fluid  theory.   In  addition  to  the  usual 

Navier-Stokes  terms,  a  finite  conductivity  term  was  included; 

obviously,  this  led  to  a  shock  thickness  dependent  on  mean 

free  path. 

Since  the  concepts  of  heat  conductivity  and  viscosity 

lose  their  meaning  at  the  high  temperatures  encountered  in 

1.  F.  de  Hoffmann  and  E.  Teller,  Physical  Review,  80, 
"Magneto-Hydrodynamic  Shocks",  692  (1950). 

K.  0.  Friedrichs  and  H.  Kranzer,  NYO-6485,  "Non-linear 
v^  Wave  Motion  in  Magneto-Hydrodynamics",  (I958). 

2.  W.  Marshall,  "The  Structure  of  Magneto-Hydrodynamic  Shock 
Waves",  Proceedings  of  the  Royal  Society,  A,  233,  357  (1956) 
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t-.hermonuclfiar  research,  a  one-fluid  theory  was  developed 
which  only  included  electrical  resistivity  as  a  dissipative 
mechanism-^  (electrical  resistivity  is  still  meaningful  for 
very  large  mean  free  path).  Although  this  is  a  sub-case  of 
Marshall's  it  is  possibly  a  better,  albeit  still  a  poor 
approximation  to  high  temperature  ionized  fluids. 

Only  a  two-fluid  theory  can  adequately  account  for 
the  phenomena  occurring  in  a  collisionles?  ionized  fluid. 
There,  if  a  shock  is  to  occur ,  the  dissipative  mechanism 
must  somehow  be  provided  by  the  cooperative  effects  arising 
between  charged  particles.   A  sufficient  condition  for  the 
presence  of  these  effects  is  that  the  Debye  length  be  small 
compared  to  any  reference  length  of  interest.   Such  a  medium 
was  named  a  "plasma"  by  Langmuir;  the  collisionless  shocks 
that  might  occur  v;ithin  it  will  be  called  here  "plasma- 
nagnetic"  in  distinction  to  the  collision-dominated  "magneto- 
hydrodynamic"  shocks. 

The  simplest  two-fluid  theories,  developed  by  Adlam  and 

h  c: 

Allen  and  Davis,  Liist  and  Schliiter-^  were  based  on  the 
assumption  of  a  zero  temperature  plasma;  they  showed  that 

5.   Rose,  M.  H. ,  "On  the  Structure  of  a  Steady  Hydromagnetic 
Shock",  NYO-7695,  (1956). 

4.  Adlam,  J.  H.  and  Allen  J.  E. ,  "The  Structure  of  Collision- 
Free  Hydromagnetic  Waves",  Phil.  Mag.  2»  ^^8-^56,  (1958). 

5.  Davis,  L.,  Lilst,  R.  ,  Schltlter,  A.,  "The  Structure  of 
Hydromagnetic  Shock  Waves:  I.   Non-Linear  Hydromagnetic 
Waves  in  a  Cold  Plasma",  Zeitschrift  fttr  Naturforschung, 
l?a,  916-36  (1958). 
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such  media  were  unable  to  sustain  a  shock  wave,  the  only 
possible  solutions  being  solitary  pulses.  The  next  simplest 
theories  for  non-zero  temperature  plasmas  due  to  Gardner, 
Grad,  et_  al_. ,  and  to  Banos  and  Vernon   ontain  essentially 
the  same  results.  More  specifically,  by  starting  from  the 
conservation  equations  with  isotropic  pressure  terms  and 
Maxwell's  equations,  they  have  shown  that  the  direction 
field  of  this  system  contains  a  closed  singularity i  thus, 
the  only  possibilities  are  a  solitary  wave  or  an  lindamped, 

periodic  oscillation.   An  interesting  result  is  that  the 

/%"  1 
wave  length  of  these  oscillations,  namely,  L  =  ]/rr^    -g 

(m  =  electron  mass,  |j,  =  permeability,  e  =  electrical  charge, 
n  =  number  density)  is  a  possible  reference  length  for  shock 
thickness. 

By  introducing  an  heuristic,  collision-dependent  fric- 
tion term  proportional  to  the  relative  velocity  of  the  two 

7 
components,  Morawetz  has  obtained  a  smooth  transition  from 

the  state  ahead  of  the  shock  to  a  damped  oscillatory  state 

behind.   As  expected,  the  shock  width  is  proportional  to  L. 

However,  since  the  real  purpose  was  to  investigate  the 

6.  C.  S.  Gardner,  H.  Grad,  et^al.,  "Hydromagnetic  Shock 
Waves  in  High  Temperature  Pla'smas"  (Papers  presented  at 
the  Second  International  Conference  on  the  Peaceful  Uses 
of  Atomic  Energy,  Geneva,  September  1958)^  NYO-2538. 

Banos,  A.,  Jr.,  and  Vernon,  R. ,  "Large  Amplitude  Waves 
in  a  Collision-Free  Plasma:  I.  Single  Pulses  v^ith 
Isotropic  Pressure".   (Physics/Plasma  Dynamics  Research/ 
LRL  -  No.  1,  Preprint  of  paper  submitted  to  II  Nuovo 
Cimento ) . 

7.  Morawetz,  C.  S.  ,  "Magnetohydrodynamic  Shock  Structiore 
Using  Friction",  NYO-8677  (1959). 
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possible  existence  of  shocks  in  infinite  mean  free  path 
plasmas,  Morav/etz  subsequently  developed  a  completely 
collision-free  theory.  By  expanding  the  particle  trajectories 
and  the  self-consistent  fields  in  terms  of  the  electron-ion 
mass  ratio,  Morawetz  has  demonstrated  to  lowest  order  the 
existence  of  a  smooth  transition  from  the  state  ahead  of 
the  "shock"  to  an  undamped,  oscillatory  state  behind. 

A  different  approach  to  the  problem  ,  described  in 
the  following  report,  uses  conservation  equations  of  mass, 
momentum  and  energy  in  which  collision-independent  stress 
terms  replace  the  usual  Navier-Stokes  terms  as  the  source 
of  a  possible  dissipative  mechanism.  This  system  of  equa- 
tions is,  therefore,  an  extension  of  the  system  with  isotropic 
stress  tensor  described  earlier  by  Gardner,  Grad,  et  al^. 
The  solutions  of  both  systems,  however,  turn  out  to  be 
essentially  the  same.   In  both  cases  the  singularity  behind 
the  shock  is  closed  which  points  to  a  pulse-like  solution. 
This  is  not  surprising  in  view  of  the  fact  that  the  equations 
are  sjnranetric  with  respect  to  the  space  variable.  The  shape 
of  the  solution  profile  is  of  interest  since  a  pulse  that 
has  a  structure  might  be  inteiTDreted  as  an  incipient  shock. 
More  precisely,  it  is  conceivable  that  the  larger  the  number 

8.   Morawetz,  C.  S. ,  "Magneto-hydrodynamic  Shock  Structure 
without  Collisions",  NYO-2885  (March  i960). 

A  preliminary  study  of  this  problem  was  described  by  the 
author  in  TID  7503  (1955). 
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of  moments  taken  into  account,  the  more  structure  v/111  appear 

in  the  pulse  with  the  result  that  the  point  where  it  returns 

to  any  initial  value  will  gradually  recede  from  the  corres- 
ponding Initial  point. 

By  using  an  expansion  method  suggested  by  C.  S.  Gardner 

it  is  shown  that,  to  second  order  in  shock  strength,  the 
pulse  has  no  structure. 
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2.   Basic  Assumptions  and  Equations. 

VJe  shall  consider  a  fully  ionized  fluid  consisting  of 
positive  ions  and  electrons  in  sufficient  nuir.bers  to  insure 
charge  neutrality  to  a  high  degree  of  approximation.   On 
the  other  hand  the  temperature  of  the  fluid  should  be  suf- 
ficiently high  for  the  collision  mean  free  path  to  be 
infinite.   Since  the  Debye  Length  is  h  =  (KT/nge2)l/2 
(K  =  Boltzmann  constant,  n^  =  electron  density,  T  =  tempera- 
ture in  degrees  Kelvin),  n^  and  T  must  also  be  of  the  right 
magnitudes  to  make  h  small  compared  to  any  dimension  of 
interest;  otherwise,  cooperative  effects  cannot  occur. 

We  assume  that  the  motion  of  the  fluid  is  confined  to 
the  xy  plane,  that  the  magnetic  field  B  only  has  a   z 
component  and  that  x  is  the  only  independent  variable  (see 
Figure  l).   This  geometry  allows  the  choice  of  a  coordinate 


Fig.  1 
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frame  moving  at  the  shock  speed  and  such  that  the  constant 
velocities  of  the  plasma  at  either  end  of  the  shock  are  In 
the  X  direction  only.   More  specifically,  the  reference 
system  is  chosen  so  that  -oo  is  "ahead"  of  the  shock  and  -t-oo 
is  "behind". 

The  motion  of  the  fluid  Is  governed  by  the  equations 
of  conservation  of  mass,  momentum  and  energy  together  with 
Maxwell's  equations.^  These  are  (using  rationalized  M.K.S. 
units). 

Conservation  of  mass: 


Conservation  of  momentum; 

J       2     XX 
(2.2)  ^  (p^Uj,   +  P^   )  =  qr(Ej,  +  Vj^) 


for   the   X  direction 


(2.^)  J.    (p    u   V     +     py^)    =   q    (E      -   u  B) 

dxrrr  r  ry  r 


for  the  y  direction 


Conservation  of  energy; 

(2.4)  ;^  (e  u  +  u  P""^  +  V  P^'""  +  q''")  =  q  (u  E   +  v  E  ) 

dx   r  r    r  r     r  r     r    ^r  r  x    r  ^• 


9.   Grad,  H.,  Notes  on  Magneto-Hydrodynamics,  IV,   NYO-64o5, 
New  York  University  (1956). 
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Maxwell's  equations 
(2.5)         curl  I  =  J 


(2.6)         curl  E  =  0 


where 


(2.7)         J=Xqr^r 


The  subscripts  "r"  may  refer  either  to  the  ions  "1"  or  to 

the  electrons  "e". 

The  notation  has  the  following  meaning: 

p_  =  mass  density  =  m  n_ 
rjp  «    r  r 

m  =  mass  of  a  particle  of  the  r   component. 

n  =  number  density  of  the  r   component. 

q  =  charge  density  (signed)  =  en  . 

u  =  X  component  of  macroscopic  flow  velocity. 
V  =  y  component  of  macroscopic  flow  velocity. 

P^"^  =  stress  tensor  =  p^5  ^  +   p^*^ . 
r  r      ^r 

p  =  pressure. 

e  =  total  energy/unit  volume  =  e  +  ^  Pr.(^T^  "^  "^r^' 


2  ^'r^^r 


e  =  internal  energy /unit  volume. 
Q^  =  heat  flow. 

E  =  X  component  of  electric  field. 
E„  =  y  component  of  electric  field. 


-  12  - 


3  =  magnetic  field. 

J  =  current. 

\i   =  permeability. 

This  system  of  equations  is  not  yet  complete  since  it 
contains  more  variables  than  equations.  The  extra  relations 
that  are  needed  are 

1)  The  perfect  gas  law  p  =  P-RTj, 

2)  Expressions  for  the  stress  and  the  heat  flow  in 
terms  of  the  fluid  variables  occurring  in  the  equations  of 
motion.   These,  therefore,  will  replace  the  usual  Navier- 
Stokes  terms.  They  are 


(?.8)  pxy  o  _l_p   i^  =  p/x 

^        '  r         2a>^  *^r  dx  r 


(2.9)  pf  -  Pr  -  sr  asr 


(2-10)  ?V  '  P.  ^  ^  '"'■ 


r         ^r        2cD     dx 

where  to     =  es/m  . 

Equations  (2.9)  and  (2.10)  lead  to  p^,  =  i-  (P^  +  P^^),  i.e., 

the  fluid  is  strictly  2-dimensional  and  therefore  has  a 

2-dimensional  stress  tensor:   P^   is  irrelevant, 

r 

The  heat  flux  vector  Q  is 

D   dRT 

(2.11)         Q^  =  (0,  qJ,  0)  where  Q^  =  |  ^  "d/  * 

r 
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since  Q^  does  not  occur  In  the  equations  of  motion,  the 
heat  flov;  term  will  be  entirely  absent.  The  derivation 
of  equations  (2.8)-(2.11)  is  described  in  Appendix  I. 

3)  The  Internal  energy  relation  for  a  2-dimensional 
gas,  namely 

With  these  relations,  our  system  of  equations  will  consist 
of  10  equations  v;ith  10  unknovms. 

Tv/o  simplifications  are  immediate;  first,  curl  E  =  0 
yields  E  =  constant  =  F.   Secondly,  the  assumption  of 
charge  neutrality  (n.  =  n  )  and  the  equation  (curl  — )x  =  J  =  0 
lead  to  u  =  u.  =  u.   Further  slnpllfi cations  are  possible 
if,  instead  of  using  separate  conservation  equations  for 
each  component,  we  use  the  total  equations  together  with 
those  for  either  one  of  the  components  (e.g.,  the  ions). 
Since  total  momentum  and  energy  are  conserved,  it  is  obvious 
that  these  equations  may  be  integrated  once. 

The  system  now  has  the  following  form: 
Mass 

(2.12)  Total     pu  =  M 

(2.13)  Ions      p.u  =  Mj_ 


-   Ik   - 


Moment  \im 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


r 


lens     ( 


Total 


^  («iV  ?lr  3^'  =  'i""-"^' 


2       p.      dv.        p        dv 
Mu  +  p  +  S-  -  ^  ^-  +  ,A-  -^  =    TT. 


■^jl  "   2co^  HF'       2cj^  ^3x" 


X 


"l^l  +  "e^  *  I  (S^  -  sf'  ak-  =   TTy 


Energy 


(2.18)         Ions 


(2.19)         Total 


d     ro  1^,22.      ^1^1   du     Pi^  ^^1 

^   [2p^u  +^Mi(u  +v^)+-_^-^33^ 


=  q^(uE  +  v^F) 


2pu  +  \  Mu2  +  I  M^v2  +  I  H^v2  +  M 

"2   ^oj.    dx       m     dx    ^      2   ^   CD.  cu^    Mx       <- 

1  e  1  e 


Maxwell 


(2.20) 


1   dB 
M-  dx 


=  qi(Vg-  v^) 


In  order  to  simplify  the  notation,  E  has  been  used  Instead 

of  E  and  B  instead  of  B,.   This  notation  will  be  adhered 
X  z 

to  throughout  the  remainder  of  the  report.  Also  p  =  p^+  p^. 
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P  =  Pa+   P  and  M  =  M.4  M  where  M^  and  M  as  vjell  as  £  * 
TT  ,  F  and  TT  are  all  constants  of  integration.  Fortunately, 
these  six  constants  can  be  reduced  to  two  by  introducing 
dimensionless  variables  and  by  using  a  coordinate  system  in 
which  TT  =  0.  The  dimensionless  variables  are 

"°^'       '1'^'      '^'rH-  ^%X'   ^l-^^l' 

-       \  -   TTx       _    TTx      _      ^± 

^e=  fr  ^e'     Pi"  — ?  Pi'   Pe"  — 2  Pe'  ^  =   > ri- 3/^  -^^ 

[vi  =  _£,       oi  =  —  =  1^     where,  since  M^  «  M.  ,  we 
e  M.  *        e  00.   M  '        e     i 

have  used  ~  =  1.     From  here  on  only  dimensionless  quantities 

^i 
will  be  used  and,  therefore,  the  bars  may  be  omitted  without 

any  possibility  of  confusion. 

In  its  final,  dimensionless  form  our  system  becomes: 


Momentum 

du  ^  ^Pi 
dx    '    dx     " 

Pi  ^s 

■    ^    dx^     ■ 

^  dp^   dv^ 
2  dx     dx 

(2-21)  ^^^--T-^-iar-asr=  ^i(^  ^  ^i^) 


(^   ^^\  "^""l  +  ^i  dfu    1  du  ^  _    /p_^) 

(2.22)         ■35r  +  ~2'—^^2^d3r-   %^^  ^^ 

dx 
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(2-2^)        "-^P+^-T^^  +  ^^=l 


(2.24)         V,  +  M^v  +J(p,--i)^=0 
^  i    e  e   2  ^^1  03^'  dx 

e 

Energy 

2 
^,,      dp.     ,,,     dv.   p. u  d  V. 


£1^1  dfu  u  ^Pi  <^^i  /l  du  ^Pi  _,  ^^  ^  ^  p^ 


dx 


2   M_v?   V?        .,     dv. 


(2.26)         2pu  +  ^  +  -|^  +  4  +  FB  -  §  (p^  _i 


Pe  ^^es  ^  1  du  ,  ^    ^ePe^   r 
e  e 


Maxwell 


(2.27)  g  =  qi(v^  -  V,) 


If  we  multiply  equation  (2.21)  by  u  and  eqi:iatlon  (2.22) 
by  v^     and  subtract  them  from  equation  (2.25)  we  are  left 
with 

^  d5r  +  2p^  S=  0   or  p^^  -^  . 

Thus  the  Ions  behave  adlabatlcally  and  the  same  Is  obviously 
true  of  the  electrons.   In  spite  of  these  adlabatlc  relations, 
the  off -diagonal  stress  terms  will  still  be  present  In  our 
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equations  because  of  the  difference  in  mass  between  the  ions 
and  electrons.   If  the  plasma  had  consisted  of  particles  of 
equal  mass,  however,  our  equations  would  have  reduced  to 
those  describing  the  isotropic  case. 


3.   Range  of  Parameters. 

We  shall  show  that  a  necessary  condition  for  the  exis- 
tence of  a  shock  is  that  the  parameters  (integration 
constants)  be  kept  within  a  certain  range  of  values.  These 
results  will  then  be  used  in  determining  the  nature  of  the 
singularities  ahead  of  and  behind  the  shock. 

Since  the  solution  we  are  seeking  must  attain  constant 
values  at  ±   co,  the  derivatives  of  the  dependent  variables 
must  approach  zero  far  from  the  shock  front.   Our  equations 
will  then  reduce  to 


(3.1) 

E  =  V.    =  v^  =   0 
1          e 

(3.2) 

F  =  uB 

(3.3)      . 

B^ 
u  +  p  +  ^  =  1 

(3.4) 

2 
2pu  +  ^  +  PB  =    £ 

Eliminating  p  between  (3.3)  and  (3- 4)  and  using  (3.2)  leads 
to 
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1 1  v/i  -  |e 

(3.5)        u  =  ^ ^— 

3/2 


hence,  it  Is  necessary  that  6  <  -r  •  Letting  e  =  ]/l  -  "f  £ 
we  obtain 

(3.6)         u=|(l±€). 

Obviously  the  "+"  sign  will  give  the  value  of  u  at  -«,  the 
"-"  sign  at  +00. 

The  lower  limit  of  £  Is  found  by  allowing  B  and  p 
to  approach  zero  ahead  of  the  shock.   It  is  permissible  to 
let  B  — >  0  since  we  are  dealing  with  fast  shocks  and  hence 
a  "complete  switch-on"  phenomenon  may  occur.    With 
B  =  p  =  0,  we  find  £  =  i  with  u  =  1  ahead  of  the  shock  and 
u  =  i  behind.  The  ranges  of  £    and  e  are,  therefore, 

I  >  €  >  0. 

We  note  that  if  £  =  §  *  ^'^^'^  ^  '^  3  everywhere,  i.e.,  we 

have  a  zero  strength  shock.   In  terms  of  our  original 

2 
variables  this  leads  to  u  =  2(2.  +  §— -)  which  states,  as 

expected,  that  a  zero  strength  shock  results  when  the 
velocity  ahead  of  the  shock  is  Just  equal  to  the  fast  dis- 
turbance speed  ahead. 

From  B  =  f/u  we  have 


10.   Cf.  report  by  Friedrichs  and  Kranzer,  pp.  30-3^>  cited 
in  reference  (1) . 
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(J.7)        B  =  |x4t 

where  the  "+"  is  to  be  used  ahead  of  the  shock,  the  "-" 
sign,  behind.  Solving  for  p  yields: 


(J. 8)  p  =  1   2  (1  ±  e).  i 1 


7^^  ''-^Tm^T^- 


since  p  >  0,  we  obtain 


Both  €  and  P  may  vary  independently  of  each  other. 
Furthermore,  each  has  a  different  meaning  since  e  deter- 
mines shock  strength  whereas  F  Is  essentially  a  measure 
of  magnetic  field  strength.   For  example,  even  with  F  =  0 
(no  B  field),  there  can  still  be  a  Jump  In  u  and  p. 
However,  for  every  value  of  €  there  Is  an  upper  limit  for 
P  which  is  determined  by  the  condition  that  p  >  0.   For 
example,  when  €  =  l/2  it  is  necessary  that  F  =  0  in  order 
to  have  p  >  0. 
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4.   Direction  Field. 

We  shall  show  that  the  system  of  equations  (2.21)  to 
(2.27)  admits  of  a  4-dlmenslonal  direction  field.  More 
specifically,  we  shall  obtain  a  system  of  four  equations 
In  v/hlch  the  first  derivatives  of  v  ,  v.  ,  u  and  B  are  each, 
separately,  expressed  In  terms  of  an  Integrated  function  of 
these  quantities. 

We  shall  now  outline  the  proceedure  for  obtaining  the 

direction  field: 

Pi  ^^1   Pe  ^^e 
1)  The  expression  -^  ^-=^  -  ^—  ^—  may  be  eliminated 

between  (2.25)  and  (2.26)  leading  to  an  Integrated  relation 

of  the  following  form: 

2 
p   -        p         p  2   M  V 

I  (Pi-  3^)b2-  p(p,-  s^)b  -(p,-  ^)  IMP,*  Pi)^  V  *  -P 

e  e         e 

V.  ^  V  p 

+  ^  +  u(l-p^-Pg-u)-  L  ]-(v^+  M^Vg)(v^p^  +  -|-H.)  =  0. 

e 

If  we  neglect  terms  In  products  of  v.  and  v  since  we  are 
dealing  with  weak  shocks  and  v.  and  v  are  both  0(€)  near 

JL  C? 

the  singularities,  this  equation  simplifies  to 

(4.1)  p^  =  i(£  +  ^  -  PB  -  u(l  -  ^)j  -  p^. 

We  retain  equation  (2.23)  which  Is  already  linear  In  dv^/dx 

and  dv  /dx. 
e 
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2)  Differentiation  of  equation  (2.24)  yields  a  relation 
for  d  u/dx   in  terms  of  first  order  derivatives.   Eliminat- 
ing d  u/dx   between  this  relation  and  equation  (2.22)  leads 
to  a  linear  relation  between  dp^/dx,  dv^/dx  and  dv^/dx. 
Since  p  ^^  l/u  ,  this  is,  in  effect,  a  second  relation 

between  dv^/dx  and  dv. /dx. 
e         1 

5)  Since  du/dx  and  ds/dx  have  the  desired  form  in 

(2.24)  and  (2.27),  we  now  have  four  linear  equations  for  the 

four  unknowns,  dv  /dx,  dv./dx,  du/dx  and  dfi/dx.  Solving, 
we  obtain 


(4.2) 


(^.3) 


p^q.         2M^  ^2 

dv,   HTpT  (^-"^)+  -T^  (1-^-Pi-  Pe-  \^ 
1  _  e'^i   

dx  J) 

2 

dv    -qi(F-^)  +  (|-  -  |)(l-^-Pi-  Pe-  I-) 

e  _  2_± 

dx  -  J) 


C*-"*)     f  =  ^I'V  "i) 


,,,     V.  +  M^v^ 
(4.5)     ^=-''    "" 


dx 


p/2 


Pe 
for  the  direction  field  where  ?  =  P^  -  ^  ^^^ 

M  p .    P  e 

D  =  I  -^ ^  (^7  -  -),  and  the  ion  pressure  p^  is  given 

by  equation  (4.1).*  In  writing  down  this  direction  field 

we  have  neglected  terms  that  are  of  O(e^),  the  reason  being 


The  electric  field  E  may  be  calculated  from  equation  (2.21); 
however,  for  o\ir  present  purposes,  we  do  not  need  to  know 
its  value. 
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that  they  will  not  be  used  at  present  in  determining  the 
nat\ire  of  the  singularities. 


5.  Nature  of  the  Singularities. 

Linearization  of  the  direction  field  in  the  neighbor- 
hood of  the  singularities  yields: 


dv.    , 

o    1         1 


(5.2)       a5r  =  r-tV(^-%)+^  (B-Bj] 

o    e         e 


(5.3)     i=^i^(ve-v,) 


(5.4) 


du     ^i+  ^e^e 


H5^ 


?y2 


where  the  subscript   "o"  may  refer  to  either  -h»  or  -w  and, 
furthermore , 

TT       _       o     o        ,   e    /      3    .    C  Ox 

i           e^i^  o           '^       u^       u^ 

o  00 


B        =  — °— °  u ^B 

Vi  -  a)gP^  o       Pq        o 

o 
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£  3  '^- 


e    ^1„    o        o 


The  equation  satisfied  by  the  elgen-values  Is 


1.     1.   e    e     oei^2 


(5.5)     X  +     q^  D^  +  p  p  /^  +  p^  d;^ D^"  /  ^  "^ 

2q. 


"o  o     o  o      o'    o      o 


.  _i0(i  ,  M^)(„^^B,^  .  „^^B,J  =  0. 


o  o 


2  1 

Letting  A  =  y  and  allowing  terms  multiplied  by  M  =  -— 

e   cOg 

to  approach  zero  (this  is  consistent  with  the  assumption 
that  the  electrons  are  adiabatic),  we  obtain: 


•  16  qj.   MgU^  /  ^    ^ 
(5.6)     y  =  |t|Va^.   ^g°,     (^|-| 

or  y  =  a  .  b  where  a  .  -^q,^  D^  ^  D^;  +  -T^^ 5^ 

and  b  =  i  /     .  Behind  the  shock,  (f/u^  -  5/2  >  0,  hence 
b  is  real  and  also  b  >  a,  (note  that  a  >  0).  Therefore, 
the  only  stable  solutions  will  be 
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lA,  =  i  )/a-E 


lAg  =  -1  /a^ 


A-,  =  -  ya+b       (with  A^,  =  /a+b  ,  unstable) 

where  A-,  and  Ap  correspond  to  a  closed  singularity  and  A^, 
to  a  saddle-like  singularity. 

Ahead  of  the  shock,  £/Uq  -  J>/2   <  0,  but  b  is  still 
real  since  £. /u  -  3/2  is  multiplied  by  M  .  However,  in 
this  case  we  obtain  0  <  b  <  a  and,  therefore,  the  stable 
solutions  are 


Ae  =  /a+b  ,     Ag  =  /a^ 


with 


A„  =  -  /a+b  ,    >^8  "^  '  ^^^^ 


unstable. 


Thus,  the  singularity  ahead  of  the  shock  has  the  nature  of 
a  saddle-point. 

In  order  to  ascertain  whether  the  closed  singularity 
can  open  into  a  spiral,  it  would  be  necessary  to  determine 
the  effect  of  the  higher  than  linear  terms.    We  note, 
however,  that  our  equations  are  symmetric  in  x,  i.e.,  they 
are  invariant  under  the  transformation  x  — >  -x,  E  — >  -E, 
V  — >  -v  with  B,  u,  p.  and  p  remaining  unchanged.   It 

11.   See,  for  instance,  Hurewicz,  "Lectures  on  Ordinary 
Differential  Equations",  Technology  Press  and  Wiley, 
(1958). 
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seems  improbable,  therefore,  that  the  closed  singularity 
can  ever  become  a  spiral.  I.e.,  that  the  solution  can 
consist  of  a  smooth  transition  from  a  steady  state  ahead 
of  the  shock  to  a  damped  oscillatory  state  behind.  This 
Is  confirmed  In  the  next  two  sections  In  which  a  method 
Is  presented  whereby  the  solution  profile  can  be  deter- 
mined analytically.   As  expected.  It  turns  out  to  be  a 
pulse. 


D.   Gardner  Transition  Curve  for  Small -Amplitude  Motion. 

The  first  step  towards  finding  the  transition  curve  is 
to  Introduce  a  new  set  of  variables  which  are  derived  from 
the  old  by  a  stretching  transfor^nation;  the  parameter  of 
the  transformation  is  e,  the  shock  strength.  Next,  the  new 
variables  are  expanded  in  a  power  series  in  e.  By  keeping 
only  terms  of  lowest  order,  it  is  possible  to  obtain  a  set 
of  differential  equations  which  are  non-linear  and  yet 
simple  enough  to  be  solved.  We  shall  show  that  the 
appropriate  solution  of  these  equations  is  a  pulse. 

The  following  calculations  are  made  easier  by  using  a 

slightly  different  set  of  dlmensionless  variables.  Let 

B^,  n^  and  u^  be  the  respective  values  of  B,  n  and  u  at 
o'   o      o  '^ 

X  =  -00.   Furthermore,  let 

The  author  is  most  grateful  to  Dr.  Clifford  Gardner  who 
suggested  the  method  and  calculations  described  in  this 
section. 
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B 


1/2 


U.{m^+   m^)n^ 


u 


Alf 


be  the  Alfven  speed  at  x  =  -oo.  The  new  dlraensionless 
variables  will  be  defined  by 


n 


"  =^  ' 


^i 

Va   = 


B-  =|-  ,    u-  =-H_ 


^   ^Alf    ' 


t     e 

v_  = 


e   u 


Alf 


P^  = 


(mi+  m^)n^u^^^ 


p:  = 


(m^+  m.^)n^u^i^ 


eB. 


X   = 


(m^+m^)u^^^ 


X. 


Lett-lng  a  =  mg/m^+  m^  and  omitting  the  primes,  we  have 


(6.1)     ^ 


^  p.  dv.     p^  dv^   ,   ^ 

(6.2)     ^  [u^u  +  p^+  Pe-d-a)^!  ^  ^  a  ^  ^  ^  ^  B^]  ^  0 


(6.5)     ^  [u„[(l-a)v,^  avj  *  (1-a)  ^  g  -  a  ^ 


^e  du 


dx 


■=  0 


dv 


(6.4)     u^(l-.)  gjl  +  (i-a)4  (5I  g)  =  n(u„-  B.) 
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(6.5)  i  =  n(v,-  V,) 


(6.6)  ua3fi+2Pig=o 


*^Pp  du 

(6.7)  .,  g/+?p^«H=o 


In  this  system,  the  equation  corresponding  to  (2.21)  has 
been  omitted  since  it  merely  serves  to  determine  E.  The 
energy  equations  (6.6)  and  (6.?)  have  been  v/rltten  down 
separately  for  each  component  (Instead  of  using  the  total 
momentum  equation  together  with  the  ion  energy  equation  as 
was  done  previously).   The  simplicity  of  their  form  results 
from  the  use  of  equations  (6.2)  and  (6.3). 

V/e  now  define  new  variables  u,  B,  p^^,  etc.,  by  means 
of  a  stretching  transformation  appropriate  for  small- 
amplitude  motion.   The  correct  transformation  is  partially 

indicated  by  relation  (5.6)  which  shows  that  the  wave- 

i/p 
numbers  X,  and  X^  scale  as  e^      for  small  €.  Thus  the  new 

1/2 
space  variable  4  will  be  defined  as  4  =  e  '^  x.  Now,  let 

the  initial  conditions  ahead  of  the  shock  be 


P          c"- 
^      ^      0          0 

Pi   =^   2     =        4 

1 

„          Po        ^0- 
Pe   =   2     =  — T 

1 

Uq  =  Co(l+0 
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where  c^  =  )/2^p^+  1/2J  Is  the  fast  disturbance  speed.  We 
shall  assi^me,  furthermore,  that 

U  =  C^{1+€)(1+€U) 

Pi  =  ("V-^^^  ""  ^Pl^ 

Pe  =  (^ir-)(l  ^  ^Pe) 
B  =  1  +  €B 

n  =  1  +  en 

and,  therefore, 

V^  =  C^(l  +  €)€^^^  v^ 
V«  =  C^(l  +  €)€^/^  V^. 

e    o^     '      e 

We  note  that  at  x  =  -oo  we  have  u=V4=v  =P4=p  =B= 

1    e    1    e 

n  =  0. 

Substituting  these  expansions  Into  the  system  (6.1)  to 

(6.7)  and  retaining  terms  of  the  two  lowest  orders  above 
the  constant  term,  we  obtain  from  (6.6)  and  (6.7): 

(6.8)  (1  +  eu)  ^+  2(1  +  ep^)  g=  0 


dp  ,— 

(6.9)  _     (1  +  €U)  ^  +  2(1  +  €P^)  gf  =  0- 
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This  yields 


(1+€P.  )  =(l  +  €p  )  =  - i^-^  =  1-  2€U  +  3€^  u^ 


correct  to  0(e  )  or 


(6.10)    p^  =  p^  =  -2u  +  3€U 


Using  (6.10)  with  the  Initial  conditions,  equations  (6.1) 
to  (6.5)  become: 

(6.1a)  n+u+€nu=0 


(6.2a)  u  +  B  +  e  |2c^  u  +  |  (c^-   l)u^+  ^ 


b2 


c^-   1  dv,  c?-   1  dv. 


-(l-a)c„(^^)g^  +  "»o(^8-'3rr  °(^'' 


0?-  1 


(6.3a)  (l-a)?^  +  avg+ii^  (!o^)du+^J2[(l.a)7^  +  aVg] 

.il^(!^l)(1.2u-B)||]=0(c2) 


,  r dv. 

(6.i^a)  u  +  B  +  €   j(l+n)(u+B)+  uB  +  0^(1-0)  ^r— 

o 
+   (l-a)(^)   ^j=   0(6^) 


(6.5a)         ^+  c^(v^-  v^)   +  ec^(H+l)(v^-  Vg)  =  b 
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We  now  assiome  that 


n  =  n^^  +  en^  +  en,  + 


U  =  U,    +   €Up   +    . . . 


V^    =   V.    +   €V.    +    . . 
1         ^1  ^2 


e         e^         eg 


B  =  B^  +  eBp  +   . . . 


To  lowest  order  we  obtain 
n^  =  -u^  =  -f(4) 

du. 


^o^V  ""^l^  ^^^  ""^^^ 


c?-  1  du. 


(l-a)7,^.  -V=  -  -^¥^  (^)  ^  = 


-1    ^=1 


l-2aN,°o"  ^ 


o     o 

To  get  an  equation  for  f{i),  we  substract  (6.2a)  from  (6.4a) 

and  Insert  the  values  of  the  various  quantities  In  terms  of 
f.  Using 

^       _   r  ,l-2a^,^o-  1 


^1       ^o     o      °o 


and 
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c2-  1 


we  finally  obtain: 


2 
(6.11)    ^  +  2Af  +  I  Af^ 

where  A  = 


5~TT2" 


a(l-a)  +  2a(l-a)(^^)-(l-2a)2^^^j 


The  solution  of  equation  (6.11)  Is  easily  seen  to  be 


f  =  -2  sech"^  -Ki        {-k   =  — 5 2.  ) 


which  represents  a  pulse. 
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7.  Second  Order  Approximation  to  the  Gardner  Transition  Curve. 

The  calculations  described  In  this  section  will  show  that 

the  second  order  terras  In  the  system  (6.1)  to  (6.7)  do  not 

introduce  any  essential  modifications  of  the  solution  found 

In  the  preceedlng  section.  More  specifically,  these  terms 

do  not  Introduce  a  structure  Into  the  pulse  which  remains 

simple. 

With  second  order  terms  Included  and  using  p^^  =  p  = 
o 

-2u  +  5€U   and  a  =  0,  eqxiatlons  (6.2)  to  (6.5)  become 
(7.1)     5  +  5  +  .  J|(c2-  l)u2  +  Sc^S  +  f  -  c^(Vi)^  j 


(7.2)     V,  .  i^(4^)||  *   c  [21^   .i-(!°:i)(X-B)^ 

+  [^  (^)B(B-l)  -  ^  (^)(1-B)]  f  1  =  0 


* 


It  had  been  previously  reported  that  these  terms  did 
Introduce  a  structure  (cf.  Bulletin  of  the  American 
Physical  Society,  Series  II,  Vol.  5,  No.  4). 
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_  _    r dv.    c'-l 

(7.3)     u+B  +  €  |(l+n)(u+B)+uB  +  c^  ^  +  (-^) 


^^1  ,  /o-\   d^u 


dv^    c?-l 


+  ^'|2c^  ^  +  (-Y)^  [(-2u-B+l)^]  +E(u..B) 
+  (n+l)\:S  J  "  ^ 

(7.4)     §1  -  c^(7^-  7^)-  c^.(K+l)(V  v^)-  €2Kc^(7^-7^)  =  0 

Equations  (7.2)  and  (7.3)  yield,  respectively: 

c^  1 
Vi  =^  (-°^)^n.e[^-.--^  ^^ 

and 


-  (-Y)  |l+€[l+B+2u]j^ 


—   —    —2    2—3 
B  =-u  +  €U   -  €  u-^. 


Substituting  these  expressions  Into  equations  (7.1) 
leads  to  the  following  second  order  differential  equation; 

(7.5)     (%)V-...]  0  -  {%)'(f  f  ^  -0^  ^ 


+  I  C^  U^  +  €(C^  u  -  2u^)  =  0 


Setting  u  =  u^"'")  +  eu^^^  +  €^u^^^  +  ...,   we  find 
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(7.6)     ^\   /^f  (2^ju(l))u(2)  ,  2,(1)'.  ,2  ,(1) 

(1)                            2                        ^^% 
where,  obviously,  u^  ^  =  -2  sech  X^   (A  =  — 7^ ). 

(=0-  1) 

We  shall  first  find  the  solutions  of  the  homogeneous  equa- 
tion. Using  the  standard  technique  for  reducing  a  second 
order  linear  equation  to  one  of  first  order,  we  set 
u(2)  /-u^^)  ^  yj(^)  and  obtain  the  Riccatl  eqviation 


(7.7)     ^  +  w2  +  k2(2+5u(=^))  =  0 
where 


HT 


We  shall  be  interested  in  solutions  valid  near  ^  =  0  and 

since  we  want  du^  Vd^  =  0  at  ^  =  0,  W   will  certainly 

2 
be  small  in  this  region.  Thus,  neglecting  W  ,  we  obtain 

an  one  approximate  solution  of  the  homogeneous  equation: 

.2^2        r..2/.2 


u(2)  =e-K^  (coshH)6K-/^%v. 


The  second  solution  w  may  be  found  by  setting  w  =  v  /  Zd^ 
where  Z(^)  is  an  unknown  function  whose  form  is  determined 
by  substituting  w  into  the  original  equation.   We  find 
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w  = 


U  = 


V  /  dz/v  . 

The  particular  solution  that  we  are  seeking  has  the  form 

c^{i)v{i)  +   C2(4)w(4)  where 

Ci(4)  =  -  /  f(4)w(Ud^ 


and 


with 


0,(4)  =  J  f(5)v(5)d5 


f(4)  =  2u 


2, ,2,      (,,,. 


'^''-  <   "<^'#)  (^ 


^^.^F"-^' 


At  4=0  it  is  easy  to  see  that  U(4)  has  a  maximum  since 
dU/d^  =  0  and  the  sign  of  d^u/d^^  is  determined  solely 
by  the  sign  of  f(0)  which  is  negative.   More  explicitly, 
we  find  that 


6k^ 


(7.8) 


dU 


~^.    GY.^ 


-2K'^4(cosh  y^i)    ^  +   ^(cosh  A-)''  sinh  A^e 


-^h^ 


.j  f(2)v(z)dz  |-|^+J^+  e-K^^^|-2K^4(cosh  A4) 
o  o  ^  ^-^^  [_  |_ 


6k' 


6k' 


-1 


o  p 

+  ~-(cosh  A4)^   sinh  A^ 


I  d4 


J7 


(O 


6k^ 
2  2        ~~? 
f(Oe"^  ^  (coshAO^  di; 
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If  Y?«   1  and  A  «  1,  both  of  which  are  satisfied  if 

2 
c  -1  »  c   or  2p  »  c  ,  then  the  various  functions  appearing 

in  (7' 8)  will  be  well  approximated  by  the  first  few  terras  in 

their  power  series  expansions.  Retaining  terms  up  to  and 

.3  * 

including  t,-^ ,   we  obtain  finally 

which  can  only  be  zero  at  |  =  0.  This  is  the  case  that 
is  of  interest  in  thermonuclear  research. 

On  the  other  hand,  if  K  »  1  and  A  »  1,  the  first 
few  terms  of  the  expansion  are  no  longer  sufficient  to  deter- 
mine the  zeros  of  du/d^.  However,  K  »  1  means  2p  «  c 
which  corresponds  to  the  cold  plasma  case  and  here  we  know 
that  the  solution  is  a  simple  pulse. 

It  is  not  yet  clear  whether  these  results  may  be 
extrapolated  to  the  intermediate  case  of  a  cool  (not  cold) 
plasma. 


The  W  term  in  the  Riccati  equation  contributes  terms  in  ^-^ 
to  du/d^;  therefore,  to  this  approximation,  w?  may  neglect 
it. 
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Appendix  .    Evaluation  of  the  Stress  Terms. 

The  equation  of  transfer,  obtained  from  the  Boltzmann 
equation  by  taking  moments  over  velocity  space,  form  an 
undetermined  system.  This  situation  may  be  corrected  by 
assiomlng  that  certain  relationships  exist  between  the 
higher  order  and  the  lower  order  moments  and  their  deriva- 
tives. The  form  of  these  relationships  will  obviously 
depend  on  the  properties  of  the  medium.   In  a  collision- 
dominated  gas,  for  example,  where  the  departure  of  the 
distribution  function  from  the  Maxwelllan  is  small,  the 
Navler-Stokes  relations  are  appropriate.   In  a  colllslonless 
plasma,  on  the  other  hand,  we  shall  assume,  in  accordance 
with  guiding-center  theory,  that  the  electric  and  magnetic 
fields  are  dominant. 

In  the  ensuing  calculations,  the  conservation  equations 

12 
will  be  used  in  the  form  derived  by  Kolodner.    We  shall 

only  need  the  following  two  equations : 

(I.l)     De  +I-.Q  +p   :|H  +  ^F   .U=o 
vj-.x;     i^c^   ax   ^rraxpr    r 


(^•2)      4^r]  ^[h   4^^}  ^  ^[P,-  1^]  -  2T,{Pr^B}  + 


P 


r 


+  ^TVr'^r]'° 


12.   Kolodner,  Ignace  I,,  "Moment  Description  of  Gas 
Mixtiires  -  I",  NYO-7980,  pg,  19  (1957)- 
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where 


U  =  u  -  u   (here  u  is  a  vector:  in  terms  of  our 
r    r      ^  ' 


previous  notation,  u  =  u,  u''  =  v) 


^r  =  f-  ^  •  ^  +(q-Yr)E  +(J-Tr.  pu)K  B] 
Fp^  =  m^n^  c  ^  c  2  c  5 


therefore 


r    r     r 


Now  assTjme  that 


E  -^OCe"-"-)    ,  B  ^  0(€"-'-) 


''      o^^^)'        ^-^°(^°) 


TT;^^''^  0(6),      TT/^--o(€),      Fp)^o(€) 

where  TTp^"^  =  P^"^^-  Pj,  and  TTj,^^  =  Pj,*^"^"  ^r'  ^^^  e  is  a 
parameter  of  sraallness.  Equating  terms  of  0(e~  ),  we  obtain, 
for  out  geometry, 

r      r     ^r 
and 

^^^  =  ^J^   =  0- 
r      r 
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Next,  equating  terms  of  0(€  )  leads  to 


^K]^^h'  m-  ''^rfr^  ^\  =  ^- 


It  Is  easily  shown  that  the  first  term  of  this  equation 
vanishes  and  that  the  remainder  yields 

.  .    ,     du     .. 

P  IJ  =  J^  n  -^  =  pJi 
r     2cD_  ^r  dx     r 
r 

P  ii  =      ^  ^ 
r     ^r   20)^  dx 

.  .        p   dv 
^r   "  Pr  ^  20)^  dx 

Since  the  relations  for  Q   may  be  obtained  in  an 
entirely  similar  fashion  from  the  transfer  equations  for 
the  next  order  moments,  we  shall  not  repeat  the  argument 
here. 
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